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Abstract 

The thermal properties of cold dense nuclear matter are investigated with chiral perturbation 
theory. The evolution curves for the baryon number density, baryon number susceptibility, pressure 
and the equation of state are obtained. The chiral condensate is calculated and our result shows 
that when the baryon chemical potential goes beyond 1150 MeV, the absolute value of the quark 
condensate decreases rapidly, which indicates a tendency of chiral restoration. 

PACS numbers: 24.10.Pa, 12.39.Fe, ll.30.Rd, ll.30.Qc, 21.65.+f 



1 



I. INTRODUCTION 



The phase transition of strongly interacting matter is of great interests to physicists. 



The expected phase diagram of neutral strongly interacting matter given by Ref. 



l| in the 



framework of N JL model is shown in Fig. [H As the temperature and the chemical potential 
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FIG. 1: The phase diagram of QCD for neutral strongly interacting matter given by Ref. 



increase, there is a transition from hadron matter to quark gluon plasma (QGP). For small 
chemical potential, experiments at RHIC have found evidences for the existence of ideal- 
liquid-like QGP whose properties__need further investigations , and LHC gives 



more information about QGP 



71]. Theoretically, lattice QCD predicts that the transition 



of hadron phase to QGP is a crossover with a critical temperature in the range 150 MeV 



lOj . The mechanism of crossover has been studied in many models 



90 MeV 151, 181, 

111 . 1121 . [131, [1^. For the low temperature and high chemical potential regime, our knowledge 



is quite limited. Experimentally, it is impossible to achieve such a condition in laboratories 
on earth. A natural laboratory holding such cold, highly compressed matter is neutron 
star about which more observations and theoretical work still need to be done 15, [l6 |. 
Theoretically, lattice QCD has difficulties in low temperature and high chemic al p otential 
regime. Model analyses predict that the phase transition is of first order QflQQ. As 
shown in Fig. [H for strongly interacting matter with low chemical potential it is in hadron 
phase and chiral symmetry is spontaneously broken. With the increase of chemical potential. 
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chiral symmetry is restored and the hadron matter transits into QGP. When the chemical 
potential increases to a certain point, QGP is expected to change into a superconductivity 



phase 19|, |20[ , the existence of which still needs to be verified. 



To understand the thermal evolution of cold dense strongly-interacting matter, one 
should first study how the baryon number density and the pressure change with the chemical 
potential and investigate the spontaneous breaking of chiral symmetry and its restoration 
where the key order parameter is the chiral condensate. However, at present an analysis 
of the thermal properties of cold dense strongly-interacting matter from first principles of 
QCD is not possible, so one has to resort to various nonperturbative QCD approaches and 



models. Chiral perturbation theory (xPT) 2l|, |22|, |23|, |2J| (for recent overviews, see Ref. 



261]) provides us an efficient method to investigate the problem in the regime where 



chiral symmetry is spontaneously broken. 

In this paper, the thermal evolution of nuclear matter at zero temperature is studied 
with xPT. In section II, the chemical potential corresponding to the baryon number density 
of QCD is introduced and is included into the effective Lagrangian of nucleon as an external 
field. This method is analogous to the one used by Gasser and Leutwyler in Ref. 22, l23 |. 
In section III, the baryon number density, susceptibility and pressure are calculated and the 
equation of state is obtained. In section IV, the thermal evolution of quark condensate is 

n 

considered. The evolution at different temperature has been studied within xPT in Ref. |27J . 
Our result at zero temperature is consistent with the low temperature result obtained there, 
and shows the trend of chiral restoration with the increase of chemical potential and gives 
a restriction on the range of applicability of xPT results. Section V gives the conclusions. 



II. CHIRAL PERTURBATION THEORY AT ZERO TEMPERATURE AND DEF- 
INITE BARYON NUMBER DENSITY 

In QCD, the baryon number density corresponding to the conserved charge of U{1)b 
symmetry is ^q~^q, where q{x) denotes the quark field operator. The baryon number of the 
system is given by {■^q^q) (the expectation value of the baryon number density operator on 
the vacuum state at definite /i). At present, it is impossible to calculate {^q^q) directly 
from the QCD Lagrangian. However, based on Weinberg's idea [2.1| that a hadron system 
could be described by an effective Lagrangian including all the possible terms consistent 
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with the assumed symmetry principles, one can calculate the vacuum expectation value of 
quark operators effectively in xPT. 

The calculation of {qjij,q), (^757/^5'), {qq) and (qj^q) was carried out by Gasser and 
Leutwyler {22, 



231] as follows. Introduce four external fields a^, v^, s and p: 

J^IT = C^'''' + qriv, + j,a,)q - q^is - il,p)q, (1) 



where Cq'"^ is the QCD Lagrangian in the chiral limit. The transformation properties of 
external fields are determined through the invariance of Cf^^ under local chiral transfor- 
mation. In the end of the calculation, the scalar external field s is set to be the mass 



matrix of quarks. According to Refs. 3, |23|] , if in the low energy regime the effective 



Lagrangian consists of all the possible terms, the results calculated from this Lagrangian 
are equivalent to those calculated from QCD. So, the generating functional deduced from 
this effective Lagrangian is equivalent to the one deduced from QCD. The generating func- 
tionals are functionals of external fields which are included in the effective Lagrangian in a 
systematic way to ensure the local chiral symmetry. Now the vacuum expectation value of 
quark operators q^^q, qi^l^q, qq, qi^q can be obtained from the derivative of the generating 
functional of the hadron system with respect to the external fields corresponding to these 
operators. 

To calculate the expectation value of the baryon number density {\q'^q) for cold dense 
nuclear matter, we use the same method and introduce the chemical potential as follows. 
Under local the quark field transforms as 

q ^ e'^'(")g, (2) 

where B is the baryon number operator. To ensure the local U{1) b symmetry of the system, 
we must introduce an external field which transforms as 

b^-^h^ + d^e{x). (3) 

Consequently the derivative acting on the quark field must be replaced by the canonical 
derivative 

d^^q d^q - B%^q = d^q - -h^q. (4) 

Set 

bf, = Sf.oHb (5) 



in the end of the calculation, then /i^ is the baryon chemical potential. Here and in the 
following, in all calculations we assume flavor SU(2) symmetry, which means equal chemical 
potential and mass for u and d quarks. The density of the conserved charge corresponding 
to U{1)b symmetry is the baryon number density operator 1^70^ = ^Q^Q- Its vacuum 
expectation value can be calculated from the derivative of generating functional 

In our approach, it is equivalent to replacing bo by fif, and then performing the partial 
derivative with respect to fib, and the result turns out to be 

dZ[fib] 

On the other hand, the nucleon field transforms as 



(7) 



tp e^^^(^V (8) 
under local U{1)b and the derivative of the nucleon field must be replaced by 

df.ip d^il) - iB^b^ip = d^tp - ibo^iibtp. (9) 
Accordingly, the pion-nucleon effective Lagrangian /^^^r for cold dense nuclear matter differs 



from the usual one in xPT |25l . [26] in that the canonical derivative of the nucleon field has 
an additional term —idofifJ-bi^- The Lagrangian for the pure pion fields Cj^ is the same as the 
one in the usual xPT, where the pions are taken as Goldstone bosons corresponding to the 
spontaneous breaking of chiral symmetry. 

According to power counting, when calculating to the order of 0{p'^), we need only to 
consider the first 4 orders of the effective Lagrangian. For pure pionic Lagrangian, only the 
leading order term C^'^ is needed: 

£(2) ^ ^Trid,Ud^U+) + ^TrixU+ + Ux^), (10) 

where U denotes the Goldstone boson fields 

U = expC^) (11) 

with 

n = I : ^"7 1 (12) 




and f being the three Pauh matrices, x = 2i?o(s + ip) with s and p being the external fields 
introduced in Eq. ([1]). In our calculation, p = and s equals the mass matrix of u, d quarks. 
Under fiavor SU{2)f, x can be expressed as 

X = 2BoM, = 2Bom = Ml (13) 

For n — N lagrangian, the terms contributing to our calculation are 

C^;l = c,{x+W - ^{u,u,){ijD>^D'^ij + h.c.) + f mX^^ + • • • , (14) 

where 

= '^[u+,d,u], (15) 

In Eq. ( IT^ . m denotes the nucleon mass in the chiral limit, Ci,C2,C3 and Ci are the low 
energy constants. is the pion decay constants. The values for the constants used in our 
calculation are shown in Table [H Here d^s and e'2S are not shown because they are irrelevant 
to our calculation. 

TABLE I: The numerical values for m, F-k^ Qa, ci,C2,C3 and ei used in our calculation. ci,C2 
and C3 are the same as the ones used in Ref. m and ei are fitted through the nucleon mass. 



M^[MeV] 


m[MeV] 


F^MeV] 


gA 


ci[GeV-^] 


C2[GeV~^] 


C3[GeV-^] 


ei[Gey-2] 


137 


896 


92.4 


1.27 


-0.90 X 10-^ 


3.3 X 10^3 


-4.7 X 10^3 


1 X 10"^ 



It can be seen that there are second and even higher order derivative terms of nucleon 

(2) 

field in the effective Lagrangian, such as the second term in C)^^. Due to the presence of 
those terms, the density of the conserved charge of U{1)b symmetry in nuclear system is 
not ip'^i/j. This is the reason why we introduce the baryon chemical potential at the quark 
level instead of adding a Ubi^^^ term directly to the effective Lagrangian. Besides, one 
may notice that the Feynman rules are different from those of the usual xPT owing to the 
additional term in the canonical derivative of the nucleon fields. 
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III. THE THERMAL PROPERTIES OF COLD DENSE NUCLEAR MATTER 



In order to study the thermal properties of cold dense nuclear matter, we first calculate 
the baryon number density and its susceptibility. According to the analysis in the above 
section, the baryon number density can be calculated in xPT from the effective Lagrangian 
of the nuclear system. The pion-nucleon effective Lagrangian can always be written as 

C^N = i^kij, (16) 

where according to Eq. (fT4l) . K is of the form 

K = i'jf.d^ + 7o/i6 -m + O (17) 

with O denoting all the remaining terms of K. The baryon number density of the nucleon 
can be calculated to be 

= DU{DetK)e'I'^'''^^ (18) 
= j DUTr{K-^^^K)e' S ^'""^^ 

where Ceff = C,-kN + i^-w- K and O denote the results of the Feynman path integral after 
integrating out the nucleon field. The most important contribution to {\q^q) from K is 

TriK~^^,). (19) 

The remaining term 

Tr{K-'^0) (20) 

comes from higher order terms in the effective Lagrangian which include the canonical 
derivative of the nucleon field, such as the second term in d^^}^. The contributions of such 
terms are far smaller than the leading term due to the smallness of C2- Therefore, we can 
neglect them in our later calculation and set 

n(/zb) = j DUTriK-^-fo)e'I'^'''^-. (21) 

After performing the Feynman path integral on the pion fields U, becomes the propa- 
gator S of the nucleon field. Then the baryon number density is 



^if^b) = I dp'Tri^oS) (22) 



It should be noted that the nucleoli propagator is different from the usual one because 
in the Feynman diagram the momentum of nucleon is always added with a term 6^Qfib- If 
the calculation is restricted to order Olp"^), i.e. we only consider the diagrams shown in 



Fig. 2, then the self-energy is the same as what was shown in Ref. [2J] except that the 
momentum p of nucleon field is replaced by p' = {po + ^biP}-, wherever it appears. The 
explicit expression for the self-energy is 



-4ciM2 + E„ + Sfc + Se + eiM^ + 0{p'' 



^ TV 



[m— 
■ C3}, 



(23) 



where d=4, / + 7o;Ufe, 

J 1 aVi-n^ 



■ArcCos{j 



U + OL 



J(l) 



2p' 
Ml 



^{b''-m2 + M2)/ + AJ, 



(24) 



87r2 



and 



(25) 



In fact, Sa^b^c coincide with the corresponding one in Ref. [2^, except that the nucleon 
momentum p is replaced by p' . The renormalization scheme used here is the one introduced 
by Leutwyler in Ref. 2J| based on the "infrared regularization" , which we denote as IR. 
The baryon number density can now be calculated to be 



n(/ib 



d'^pTr{ 



7o 



(26) 



(27r4) 7 'j6'-m-S(i)') 

Remember that the chemical potential is introduced in Euclidean space, po is imaginary 
and the integral is of the form 



d'^p = / d^p 



dpo. 



(27) 



Besides, note that the baryon number density vanishes at zero chemical potential, so if we 
let Po = Po + then 



(2^4 



dV( / dp'o - 



dp'o)Tr{ 



7o 



^' — m — S(p') 



(28) 
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FIG. 2: The one-loop diagrams contributing to the self-energy of the nucleon. The crossed vertices 



denote contribution from £ 
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ttN 



lm[pj i 




1 


\ 






1 


\^ 


Re[pJ 



FIG. 3: The integration path of po for the baryon number density of the nucleon in Eq. (|28 



The integration path in the complex plane would be closed as was shown in Fig. [3l According 
to Cauchy integral formula, the value of the integral is determined by the analytical prop- 
erties of the integrand in the region encircled by the integration path. However, a direct 
application of the IR results 2J] in the full complex po plane is illegal. In IR, the loop 
integral H is divided into two parts, the singular part I and the regular part R, and the 
R part is dropped out in renormalization because it does not fit the infrared singularity of 
the integral. Considering that the analysis of the singularity is based on an expansion of 
the nucleon momentum around the nucleon mass 171^-, when the nucleon momentum goes 
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far beyond niN-, the expansion is illegal and so the regular part R can not be dropped out. 
In the regime where the momentum of the nucleon is much larger or smaller than niNi the 
result of the loop integral should he H = I + R, not I. 

According to above consideration, the integrand in Eq. (1251) has three poles and one 
branch cut in the complex Pq plane with the real part of po positive and not larger than 
1200 MeV. The regime beyond 1200 MeV is out of the consideration of xPT because there 
the spontaneous breaking of chiral symmetry is expected to be restored. Then the poles 
and the branch cut relevant to our calculation are shown in Fig. |H Obviously, the pole 



♦ pole I 



Re[pJ 



lm[pj 



FIG. 4: The poles and the branch cut of the integrand in Eq. (|28p in the complex po plane. 
The positions of the three poles are at po = ^p^ + p^, where n denotes a, b and c. pn for 
each pole is pa = 938, ph = 1152 + i337, pc = 1152 — i337. The branch cut starts at d where 
PQ = a/f^ + {m + M^f. 

denoted as "a" on the real axis corresponds to the ground state of the nucleon (N938), and 
the conjugated pairs "b" and "c" can be explained as the excited modes of the nucleon. The 
brach cut starts at po = \J'P^ + (m + M^)^. 

The results for the baryon number density are shown in Fig. El where the contribu- 
tions from poles and branch cut and the total contribution are shown respectively. The 
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FIG. 5: The curves for the baryon number density of cold dense nuclear matter. The left top shows 
the contributions from poles. The right top shows the contribution from the branch cut. The left 
bottom shows the total density and the right bottom shows the total density in the unit of the 
saturated nuclear density po- 

susceptibility of baryon number density x is calculated to be 

= (29) 

The pressure P of the nuclear matter can be calculated through the relation 

dP 

n{l^,) = (30) 

which comes into existence when the temperature of the system stays the same. The results 
of X ci^iid P are shown in Fig. O The P — n curve for the equation of state is shown in Fig. [71 



The results indicate that the baryon number density n(/i) and pressure P of the nucleon 
are zero for baryon chemical potential /Ufe smaller than the mass of the nucleon 938 MeV. 

11 




FIG. 6: The curves for the pressure and the baryon number susceptibihty of the cold dense nuclear 
matter. The left shows the baryon number susceptibility. The right shows the pressure. 




FIG. 7: The curve for the equation of state of cold dense nuclear matter. 

That is, /i = 938 MeV is a singularity. This result agrees qualitatively with the general 
conclusion of Ref. ^]. In that reference, based on a universal argument, it is pointed 
out that the existence of some singularity at the point fi = fiQ and T = is a robust and 
model-independent prediction. A recent model calculation using the rainbow approximation 
of the Dyson- Schwinger approach [29] also support this. When is larger than 938 MeV, 
77, (yu) and P increase with the increase of baryon chemical potential. And when jif, is larger 
than 1152 MeV, the rate for the increase of and P is slowed down. The peak in the 
curve of susceptibility indicates a second order phase transition around n{fi) = 1150 MeV, 
which is due to the exited modes of nucleon. Physically, the critical point for ri(yu) and P to 
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go from zero to a nonzero value should be smaller than 938 MeV due to contribution from 
the binding energy. However, our calculation cannot reveal such contribution because the 
pion loops are not considered. It can be expected that when calculating to higher order, the 
results would be better and the quenching effect would be revealed. 

Remember that xPT is available only when chiral symmetry is spontaneously broken. 
As we have indicated in the introduction, with the increase of chemical potential there is 
a transition from hadron phase to QGP for which the chiral symmetry is restored. So, the 
applicability of above results should be limited in a certain range. In order to analyze the 
restoration of chiral symmetry and determine the range of applicability of our results, we 
need to calculate the quark condensate {qq) which characterizes the spontaneous breaking 
of chiral symmetry. This will be done in the next section. 



IV. THE EVOLUTION OF THE QUARK CONDENSATE 

According to C^^t^ in Eq. ([T]), the quark condensate {qq) can be calculated from the 
derivative of the generating functional with respect to the corresponding external fields 

m = Q-J D,pDijDUe'I''"^^^ffl=M,. (31) 

As the baryon chemical potential is not included in £7^, the pure pion loop does not contribute 
to our calculation. By analogy with the deduction in the last section, the result comes out 
to be 

where rh denotes the mass of u and d quarks under flavor SU(2) symmetry. The integral is 
performed in the same way as in Eq. (!28ll . 

d'v = / d'p{ / dp',- dp',) (33) 

J J —ioo+fii, J— too 

with the integration path shown in Fig. [3l However, since the quark condensate does not 
vanish at zero chemical potential, the result of the integral in Eq. (l32l) is not (qq) but instead 
{qq) — {qq)o, where {qq)o denotes the quark condensate at zero chemical potential. 

The curve for —{{qq) — {qq)o) is shown in Fig. [8] and the curve for the quark condensate 
— {qq) is shown in Fig. [9l {qq) remains the same as the quark condensate at zero baryon 
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FIG. 8: The curve for the evolution of —{{qq) — {qq)o)- 
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FIG. 9: The curve for the evolution of —{qq) 



chemical potential, i.e., (qq) = {qq)o- A recent model calculation using the rainbow approx- 
imation of the Dyson-Schwinger approach gives the same result 30|. For /i^ larger than 
rriN, the absolute value of the quark condensate tends to decrease as //f, increases, and the 
decrease speeds up when /i^ is larger than 1150 MeV, which means that the spontaneous 
breaking of chiral symmetry tends to get restored. Therefore, the xPT result is applicable 
only for fib smaller than 1150 MeV where the spontaneous breaking of chiral symmetry still 
works. Our results are consistent with the HB^PT result at low temperature in 2?!, where 
the evolution of the quark condensate is investigated for different temperatures. And in 
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principle our results are also consistent with that in of Refs. [3l|, |32]. Although xPT is not 
applicable, the point at which the hadron-QGP transition occurs can be estimated to be 
around 1150 MeV, which is consistent with the result of //g = ~ 365 MeV expected in 



NJL model 



V. CONCLUSION 

The thermal properties of cold dense nuclear matter are analyzed and the results for 
the evolution curve of the baryon number density, the baryon number susceptibility and the 
pressure are shown in Fig. [5l [61 The curve for the equation of state is shown in Fig. [71 
An analysis of the quark condensate shows that chiral symmetry tends to get restored 
for chemical potential larger than 1150 MeV. Our results are applicable only in the range 
fib < 1150 MeV, where the spontaneous breaking of chiral symmetry approximately holds. 

Our results show that the transition of nuclear matter to QGP occurs at around 1150 MeV, 

ft 

which is consistent with the prediction in [li|. To understand the full thermal properties of 
cold dense strongly interacting matter, we need to investigate the region above the phase 



33, 



34 



351 ]. which is apphcable 



transition point using a method, such as the DSE method in 
for cold dense quark matter. As is indicated in section III, the contribution of binding energy 
and the quench effect are not included in our calculation. To improve our result we need to 
make the calculation to order of 0{p^). 
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